Family of Solvable Generalized Random Matrix Ensembles with Unitary Symmetry 
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We construct a very general family of characteristic functions describing Random Matrix Ensem- 
bles (RME) having a global unitary invariance, and containing an arbitrary, one- variable probability 
measure which we characterize by a 'spread function'. Various choices of the spread function lead to 
a variety of possible generalized RMEs, which show deviations from the well-known Gaussian RME 
originally proposed by Wigner. We obtain the correlation functions of such generalized ensembles 
exactly, and show examples of how particular choices of the spread function can describe ensembles 
with arbitrary eigenvalue densities as well as critical ensembles with multifractality. 
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The concept of Random Matrix Ensembles (RME), 
originally proposed by Wigner to describe statistical 
properties of the eigenvalues and eigenfunctions of com- 
plex nuclei Q , have proved to be a very useful idea in the 
studies of a wide variety of physical systems including 
equilibrium and transport properties of disordered quan- 
tum systems, quantum chaos, two-dimensional quantum 
gravity, conformal field theory and chiral phase transi- 
tions in quantum chromodynamics as well as financial 
correlations and wireless communications 0, 1^ ■ The un- 
derlying reason for such a wide range of applications is 
the universality of the correlations between eigenvalues 
of given classes of RMEs. For example, once appropriate 
variables are chosen in which the mean spacing between 
eigenvalues is unity, the Gaussian ensembles of random 
N X N matrices (with given symmetries) have, in the 
N ~* oo limit, a universal (zero parameter) form for 
the nearest neighbor spacing distribution or the spectral 
rigidity (number variance of eigenvalues in a given range) 
known generally as the Wigner-Dyson distribution 
The appropriately scaled energy levels of complex nuclei 
and transmission levels of weakly disordered mesoscopic 
metals both follow the above universal distributions even 
though the physical sizes of the systems differ by about 
nine orders of magnitude j3| • 

More recently, much interest has been generated in 
finding RMEs that deviate from the universal properties 
of the Gaussian Ensembles in specific ways. One partic- 
ular example is the attempt to find 'critical' ensembles 
relevant for systems at the critical point of e.g. the Ander- 
son transition in disordered conductors where the spacing 
distribution or the spectral rigidity is known to deviate 
from those of Gaussian RMEs |5|,y,|3- Another example 
is the attempt to find a one-parameter generalization of 
the Gaussian RME that can describe a monotonic change 
from the universal Wigner-Dyson distribution to a com- 
pletely uncorrelated Poisson distribution as the parame- 
ter is changed, which may be relevant for a crossover from 
a chaotic to an integrable system [^. Yet another exam- 
ple is the observation for financial cross-correlation ma- 
trices in which statistics of most eigenvalues agree with 
the universal predictions of Gaussian RMEs but there 



are deviations for a few of the largest eigenvalues [9j. It 
is therefore of great interest to a wide variety of areas 
and disciplines to study RMEs that are in some sense 
generalizations of the Gaussian RMEs. 

In seeking generalizations of Gaussian RMEs, suit- 
able for arbitrarily many variables, it is natural to begin 
with the probability distributions Pn{X), where X de- 
notes an N X N Hermitian matrix. Gaussian ensembles 
have centered distributions that are in fact (exponential) 
functions of the single variable tr{X'^). It is natural to 
restrict our generalizations to probability distributions 
Pjy(X) = VFAr(ir(X^)). One such example that might 
come to mind could be Pjv(X) c>c exp(— [ir(X^)]^), but 
although such a proposal is satisfactory for any finite 
N, it fails to generalize to a valid new distribution as 
N oo. The clue to discovering the proper class of 
generalizations is to work not directly with the distri- 
bution Pn(X) themselves, but with their Fourier trans- 
forms, i.e., with the associated characteristic functions 
[i3 Cn{T) given by 
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(T) ^ Je' *'-(^^)Pjv(X)dI/x 



(1) 



with Cn{0) = 1. The integration is over the invariant 
Haar measure that preserves hermiticity. In the present 
work we first prove that if Cn{T) is a function of tr{T'^) 
only, then the most general Cn{T), valid for arbitrarily 
large N, can always be written as 



Cjv(r) = 



f{b)db , 



f{b)dh^l (2) 



where /(&) is any non-negative function, which may be 
chosen phenomenologically. While Eq. (0) is evidently a 
positive superposition of Gaussians, a vast family of dis- 
tributions is thereby included, e.g., Cauchy, Levy, etc.; 
this result generalizes the specialized and more restric- 
tive examples of We then show that the n-point 
correlation function for the corresponding unitary RME 
of the matrices X with eigenvalues Xi can be written 
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down exactly as 







dbfjb) 
(46)"/2 
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Here we have defined Xi = Xi/2Vb, and K^{xi, xj) = 
'^n=o 'Pri{xi)ipn{xj) is the weU-known 2-point kernel of 
the Gaussian RME llj, where (pn{x) cx e^^'/^if„(x) are 
orthonormal functions associated with the Hermite poly- 
nomials Hn{x). In particular, the one-point function is 
just the density of levels, given (in the large N limit) by 



(7Ar(.T) = 



2tt 
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dbfjb) 
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v/1 - xySNb, (4) 



where we have used the known result that for large N, 

for 



K^ix,x) 



jv(x) = ^/2N^^/^: for |a;| < V2iV and 
zero otherwise It is clear that one can obtain dif- 
ferent x and N dependence for the densities by choosing 
different f{b). Similarly, the two point cluster function 
defined as T2{xi,X2) = —R2{xi,X2) + Ri{xi)Ri{x2) has 
the form Ta = - ST2, where 



T^{XI,X2) 



4b 



t:i,X2 



(5) 



and 



ST2 = 



dbfjb) 
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o-^(a;i)CT^(a;2) - o-Ar(a;i)crAr(x2). (6) 



We call f{b) the 'spread function'. Note that the Gaus- 
sian RME corresponds to the choice f{b) = S{b — bo), 
for which 6X2 is identically zero. Other choices of the 
spread function can describe a variety of possible gen- 
eralized RMEs. Alternatively, correlation functions of 
a physically relevant RME of matrices X characterized 
by a given Cn{T) can be obtained exactly if the corre- 
sponding spread function fib) can be identified. For ex- 
ample, Cn{T) = e^^'o^*'' will correspond to a choice 
of f{b) oc 6~^/^e~''o/'', for which the exact correlation 
functions can be easily written down. 

The universal features of Gaussian RMEs arise in the 
N 00 limit and when the variables are chosen in which 
the mean level spacing is unity (this is known as 'unfold- 
ing'). In place of a;i and X2 one defines new variables p 
and ^ such that dp = a{xi)dxi^ = a{x2)dx2^ and the 
new cluster function 



Y2{p, QdpdC, = T2{xi,X2)dxidx2 



(7) 



is well defined everywhere in the limit N 00. For 
Gaussian RMEs, there exists the sum rule that the inte- 
gral of F2(p, C) over p is always unity. It has been argued 
\\'A. that the violation of this sum rule is a signature 
of critical ensembles, where the deficit of the sum rule 
(for translationally invariant cluster functions) 



rj — 1 



Y2{p,C)dp 



(8) 



is related to the multi-fractality of wave functions at the 
critical point We will show that there are choices 

for f{b) for which the sum rule is violated. These choices 
would then correspond to critical ensembles. 

We begin by proving that if C{T) is a function of 
tr (T^), then the most general Cm{T), vahd for arbitrar- 
ily large N , can be written in the form Eq. Note that 
if C{T) is a function of tr (T^) only, then from Eq. (|TJ 
we have Pn{U^ XU) = Pn{X) where U is unitary. Thus 
the distribution is invariant under a rotation of basis. 

Proof of Eq. 0): We suppose that 



C{T)^Ei\\T\\l)^ I , 



dp{X) 



(9) 



for all Hermitian T with ||T||2 = tr{T'^) < 00, where p, is 
a suitable probability measure. It follows that C{T—S) = 
E{\\T — 5112) is a real, continuous function of positive 
type, and the GNS Theorem ensures that there exist 
vectors \T) and \S) in a separable Hilbert space such that 

{S\T) ^ E{\\T ^ S\\l) . (10) 

A family of operators V{U) may be defined by 

{S\V{U)\T) = {S\T + U)= E{\\T + U-S\\l), (11) 

and it readily follows that V{U) is an Abelian group of 
unitary operators for which the operator norm || V^(C/)|| = 
1. Now consider the sequence C/(^^) = {Ul-f^} where 

f/^f ) = ^ SrAiSsM , u>0, 1 < M < 00 . (12) 

It follows that the weak operator limit 

lim (S\ViU^''y)\T)^Ei\\T-S\\l + u) ^ (S\Aiu)\T) , 

M—^oo 

(13) 

an expression that defines the operator A{u) for all u > 
0. Clearly, A{uy = A{u), A{u)A{v) = A{u + v), and 
||^(w)|| < 1, and therefore A{u) — e^"^, where = 
B >0. Hence, 



Eiu) ^ (0|e-"-«|0) = / e 
'0 



-ub 



dm{b) 



(14) 



where to is a probability measure, and the latter relation 
follows from the spectral representation for B. Finally, if 
we assume that m is absolutely continuous and replace 
u by ||r||2, we recover Eq. This completes the proof 
of Eq. I©. 

We next show that given Eq. the 71-point correla- 
tion function is given by Eq. 

Proof of Eq. 0): Given the characteristic function 
Eq. 12 the probability density is, from Eq. 

Pn{X)o, dbf{b) Je-''^^^'''>e-''<^'UVT. (15) 
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The integral over N'^ independent elements of T results 



P.Woc^°°|^e-(-^)/-. (16) 

The joint probability distribution of the eigenvalues of X 
is then given by 

P^({xJ)(x r^n(-,-x,.)^e~E-?M^ (17) 

where Jli<j(-^« ^ ^i)^ j^^^ the Jacobian of trans- 
formation from the matrix element to the eigen- 
value/eigenvector coordinates. We now take advantage 
of the known results for Gaussian Random Matrix En- 
sembles by noting that before the 6-integral, a change of 
variables Xi — Xi/2\/b changes the distribution to exactly 
Gaussian RME with some additional 6-depcndent terms; 
this leads directly to Eq. 0. This completes the proof 
of Eq. ©. 

As an example, let us consider the spread function 
/(&) oc foA^V2+--ie-£V4fce-7\ Then 



Pn{X) 



+ tr jX^) 
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tr (X2)) 



(18) 

where K^, is a modified Bessel function. In the limit 
7-^0, choices of v would include ensembles of Levy 
matrices li v = ~n where n is a positive integer, 

then in the same limit 7 — > this gives rise to Pn{X) oc 
(e^ + tr (X-^))^", while in the opposite limit 7 — > 00 we 

get Pn{X) cx {e'^+tr (^X^))-(n/2+l/A)^-^-<[e-+tr (X^))^ 

The 2-point correlation function for all 7 and s can be 
written down exactly from Eqs. (0) and 

As an explicit example of how novel behavior of the 
correlation functions may arise for some choices of f{b) 
in the large TV limit, let us consider the level density for 



fib) 



26- 



= CAT 

= 0, 



[&(5+l)]3/2' 

b>bo 



b < bo; 



(19) 



where b = y/8Nb and bo is determined from the nor- 
malization condition. We choose cn such that 60 = 
^/8Nbo > 1. Then the density obtained from Eq. 
is given by 



. ^ 1 Cn 



TT .Jx{x + 1) ' 



(20) 



which is similar to that satisfied by the transmission 
eigenvalues in disordered conductors and is known to lead 
to deviations from Gaussian RME 0, 0| . Note that the 
density diverges as l/Vx in the hmit x ^ 0. It is clear 
that by appropriately choosing the spread function /(6), 
one can obtain a variety of densities that can go to zero. 



a constant or infinity as a function of 0. One 

can also choose f{b) to obtain an A^-dependent density 
at the origin that goes to zero, a constant or infinity in 
the N ^ 00 limit. An A'^-independent finite density at 
the origin was conjectured to be important for the novel 
properties of the g-Random Matrix Ensembles reviewed 
in Q. 

In particular, f(b) can be defined such that in the limit 
A^ ^ 00 both (j{x) and ST2{x,y) become constant. In 
such cases if we define 



0(A) = / dbf{b)e 
Jo 



(21) 



then (t>{0) = 1, a{x) = cto = -^^'(0) and 5T2{x,y) = 
£^W'iO) - (</''(0))2]. We now choose p = aox, C = cjoy 
to scale the density to unity. Then the scaled cluster 
function defined in Eq. {Tj) has the simple form Y2 (p, C) = 
T2{xi,X2)/cro- Using Eqs. (jSj and © to define 1^2° and 
6Y2, we find 



1 f° 



dbfjb) . 




2b (To 



, (22) 



where we have used the known large A^ behavior for 
K^{x,y) Q, and 



SY2ip,0 



'(0) 



1 



(23) 



(0'(O))2 

The deficit of the sum rule Eq. © then takes the form 



7?=1 



N/2 



N/2 



dpY^{p,0+N 



'(0) 



(^'(0))^ 



- 1 



(24) 



If (To Of Va", the argument of the sine function in Y2 is 
independent of A^. Then the integral of Y2 is just unity 
for all f{b) in the limit N ^ 00. Using the inequality 



dbfib)b- 



-1/2 



< 



dbf{b)b-^ / dbfib) (25) 



and the normalization of f(b), we obtain ?7 > 0. The 
equahty sign holds only if /(6) cx S{b — bo), which is the 
Gaussian RME. As an explicit example of finite positive 
77 in the large A^ limit, let us choose 



fib) oc 6"e~''/^e~'^^. 



(26) 



With a = —3/4, this immediately leads to 77 = N/2^/J]j. 
If \/]3j oc Af, we get a well defined critical ensemble in 
the thermodynamic limit. Note that for finite ry, Eq. H23|l 
gives dY2 — rj/N in the N ~> 00 limit. 

The variance S](s) = (n^) — (n)^ of the number of eigen- 
values n in an interval (— s/2, s/2) is given by 



E(s) = dp da5{p - C) - Y2{p - C)]. (27) 

J-s/2 J-s/2 
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For Gaussian RME, E(s) oc Ins for large s. For critical 
ensembles discussed in the literature H, 0, S(s) cx s. 
Our choice of the spread function f{b) in Eq. which 
gives rise to a constant SY2, produces a term S(s) (x 

in addition to terms S(s) oc Ins from However, 
as mentioned earlier, this term is proportional to rj/N, 
which vanishes in the large N limit for finite 77. Eq. H26() 
therefore corresponds to a novel kind of critical ensemble 
with E(s) oc Ins for large s. 

Note from Eq. lO that in general ST2 is not transla- 
tionally invariant; clearly the critical ensembles in such 
cases will not give rise to number variance S(s) oc Ins, 
and will correspond to a different class. We also get a dif- 
ferent class of critical ensembles if the contribution to rj 
from Y2 is different from unity. For choices of /(&) which 
diverge as & ^ 0, changing the order of the integrals over 
p and b in Eq. (jS]) may not be justified, and the integral 
in Eq. (|22|l may depend on the choice of such f{b). 

In summary, we have constructed a generalized Ran- 
dom Matrix Ensemble whose characteristic function con- 
tains an arbitrary non-negative spread function /(&) with 
the only condition that f{b)db = 1. The correlation 
functions of the generalized ensembles are exactly solv- 
able for any given f{b). Various choices of f{b) lead to 
a variety of possible density of levels (Tjv(a;), which can 
depend on a; or in a variety of different ways, leading 
to possible deviations from Gaussian Ensembles. In par- 
ticular, we showed that it is possible to choose forms of 
f{b) that lead to violations of the sum rule for the scaled 
2-point cluster function Y2 {p, () where the deficit of the 
sum rule 77, as given in Eq. JHjl, is a characteristic of crit- 
ical ensembles with multifractal wave functions. Unlike 
critical ensembles discussed in the literature, these sum 
rule violations can correspond to different forms for the 
number variance S(s), corresponding to different classes 
of critical ensembles. These solvable generalized ensem- 
bles should therefore be of interest in a wide range of ar- 
eas where Random Matrix Ensembles play an important 
role. While there are only a few known examples of gener- 
alized RMEs for which correlation functions can be eval- 
uated exactly 0jl3|' ^ given model of Cn{T) or Pjv(X) 
characterizing a physically relevant generalized RME be- 
comes exactly solvable if the corresponding spread func- 
tion f{b) can be found. This opens up the possibility 
to obtain exact results for a variety of interesting and 
physically useful generalized RMEs. 
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